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The theory of linear response to perturbations of the equilibrium state, or linear 
response theory, is the subject of this series of lectures. 

Ordinary matter, if left alone, will sooner or later attain an equilibrium state. This 
equilibrium state depends on the temperature of the environment and on external 
parameters. External parameters may be the region of space within which a certain 
number of particles are confined, mechanical stress, or the strength of an external 
electric or magnetic field. 

If temperature or the external parameters change slowly enough, the system can 
attain the new equilibrium state practically instantaneously, and we speak of a 
reversible process. On the other hand, if the external parameters vary so rapidly 
that the system has no chance to adapt, it remains away from equilibrium, and we 
speak of irreversibility. 

The most important application is optics. There is a medium which is exposed to 
an electromagnetic wave. The electric field changes so rapidly that matter within 
a region of micrometer dimensions cannot react instantaneously, it responds with 
retardation. We shall work out the retarded response in linear approximation. 
There are quite a few general and important results which hold irrespective of 
a particular Hamiltonian, such as the Kramers-Kronig relations, the fluctuation- 
dissipation theorem, the second law of thermodynamics, and Onsager's relation. 
We discuss various electro- and magnetooptic effects, such as the Pockels effect, 
the Faraday effect, the Kerr effect, and the Cotton- Mouton effect. 
We also treat spatial dispersion, or optical activity and indicate how the theory is 
to be developed further in order to handle the non-linear response as well. 
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1 MAXWELL EQUATIONS 



1 Maxwell equations 

In this serious of lectures we will study the interaction of a rapidly oscillating 
electromagnetic field with matter. Therefore, a recollection of basic electrody- 
namics seems to be appropriate. 

1.1 The electromagnetic field 

The electromagnetic fields E = E(t,x) and B = B(t,x) are defined by their 
action on charged particles. The trajectory t — > x{t) of a particle with charge 
q is a solution of 

p = q(E + v X B) , (1.1) 

where v = x and p = mv / y/l — v 2 /c 2 . The electromagnetic field is to be 
evaluated at the current particle location t,x(t). 

The electromagnetic fields act on charged particles, as described by the Lorentz 
formula (1.1), and charged particles generate the electromagnetic field. dQ = 
dV p(t, x) is the amount of electric charge in a small volume element dV at x 
at time t. Likewise, dl = dA ■ j(t,x) is the charge current passing the small 
area element dA from the back to to the front side. 
Maxwell's equations read 

e V ■ E = p and —V x B - e E = j (1.2) 
Mo 

as well as 

V -B = and Vx£ + B = . (1.3) 

The first group of four equations describe the effect of electric charge and cur- 
rent, the second group of likewise four equations say that there is no magnetic 
charge (magnetic monopoles). 

It is a consequence of Maxwell's equations that charge is conserved, 

p + V-j = . (1.4) 

1.2 Potentials 

Stationary fields decouple. 

e V -E = p and Vx£ = (1.5) 
describe the electrostatic field, 

—V x B =j and V • B = (1.6) 
Mo 



1.3 Field energy 
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the magnetostatic field. 

The electrostatic field may be derived from a scalar potential, 

E = -V<p (1.7) 

which obeys Poisson's equation 

-e A(f } = p . (1.8) 

The magnetostatic field can be expressed by 

B = VxA (1.9) 

in terms of a vector potential A. Adding the gradient of an arbitrary scalar field 
A, A' = A + VA, does not change the induction field B. This ambiguity allows 
to subject the vector potential to a gauge, e.g. the Coulomb gauge V • A = 0. 
The three components of the vector potential the obey a Poisson equation each, 

-—AA = j . (1.10) 

The full set of Maxwell equations are solved by 

E = -V(f> - A and B = V x A . (1.11) 
If we now impose the Lorentz gauge 

e <i> + — V-A = , (1.12) 
Mo 

the following equations result: 

€ n<j) = p and — □ A = j . (1.13) 
Mo 

The box, or wave operator is <9q — A where do is the partial derivative with 
respect to time, divided by c which is defined as c = 1/^/eo/^o- 

1.3 Field energy 

The potential of a point charge q resting at y is the Coulomb potential 

<h(x) = ^^ . (i.i4) 

47reo \x — y\ 

If charges qi, q2, . . ■ are brought from infinity to their locations at xi, a?2, • • ■ the 
following work has to be done: 
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1 MAXWELL EQUATIONS 



For a smooth charge distribution p this expression may be rewritten into 
d 3 x"d 3 x'p(x")p(x') 1 



W 



1 f d 6 x"d 6 x'p(x")p(x') 1 f „ . , , 



By partial integration we obtain 

r.2 



w 



f*x*§-. (1.17) 



We interpret e$E 2 /2 as the energy density of the electric field. 

Using Maxwell's equation we may prove the following balance equation: 

r) + V-S = r/* (1.18) 

where 

77 = + — , S = — E x B and rf = -j ■ E . (1.19) 

2 2/jo /Jo 

As explained before, e^E 2 /2 is the energy density of the electric field. The 
energy density of the mag netic field is B 2 /2p . The energy current density 
S = (E x B) / po is also called Poynting's vector. 

rf = —j ■ E describes the production of field energy per unit time per unit 
volume. Field energy is created if charge runs counter to the electric field 
(Bremsstrahlung). Field energy vanishes if charges run with the direction of 
the electric field (Ohm's law, Joule's heat). 



1.4 Polarization and magnetization 

Q = f d 3 x p(x) is the charge of system. By p = f d 3 x x p(x) we denote its 
electric dipole moment. The dipole moment does not depend on the choice of 
the coordinate system origin if its charge vanishes. 

Denote by P the density of a probe's dipole moments, its polarization. For 
ordinary matter, which is locally neutral, this is a well defined quantity. One 
can show that — V • P is the charge density causing the polarization and that 
P contributes to the current density. 

We likewise define magnetic dipole moments m and their density M, the mag- 
netization. V x M also contributes to the current density. The charge and 
current density should therefore be split into 

p = -V ■ P + p f and j = P + VxM + f . (1.20) 

The remainders p f and jf f are the charge and current density of free, or mobile 
charges, as opposed to bound charges. 

We introduce auxiliary fields D = €qE + P (dielectric displacement) and H = 
B/po — M (magnetic field strength). They obey the following equations: 



V D = p i and V x H - D = / 



(1.21) 



1.4 Polarization and magnetization 

V B = and V x E + B = 

remain unchanged. 



8 



2 A SIMPLE MODEL 



2 A simple model 

For warming up, we develop a very simple model of the dielectric susceptibility. 
Consider an atom which is exposed to an oscillating electric field E = E(t) along 
the z-direction. A particular electron of this atom will be forced to oscillate. 



2.1 Equation of motion 

We assume that the electron's equilibrium position is z = 0. There is an equilib- 
rium restoring counter force which, for small deviations, will be proportional to 
z. We also allow for friction which we assume to be proportional to the electrons 
velocity z. Therefore, the following equation of motion has to be solved: 

z + 2rjz + VL 2 z = a(t) where a{t) = -^-E(t) . (2.1) 

m and — e are the electron's mass and charge, tu^Iq is the spring constant of 
the linear oscillator, 2mr] the friction coefficient. 



2.2 Green's function 

Because there is a linear relationship between position z and acceleration a we 
should write 

z(t) = J dsr(t,s)a(s) , (2.2) 
where r is a Green function. 

Because the coefficients of the differential equation (2.1) do not depend on time 
t one has r(t, s) = r(t — s). This implies 

z(t) = J dsr{t- s)a{s) = J dTr{T)a{t-T) . (2.3) 

r is the time difference between cause (accelerating force, or a) and effect (de- 
viation from equilibrium position, z). Causes must be earlier than their effects, 
therefore 



/t poo 
dsr{t- s)a{s)= / dTr{T)a{t-T) . 
-oo JO 



(2.4) 



r is an age, it cannot be negative. 

We differentiate (2.4) with respect to time, 



z(t) = T(0) a(t) + f dst{t-s)a{s) , (2.5) 

J — oo 



and likewise for z(i). 



2.3 Susceptibility 
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Inserting these expressions into (2.1) results in 

f + 2rir + nlr = o (2.6) 

and 



r(0) = and r(0) = 1 . (2.7) 

This differential equation with constant coefficients is solved by r oc e ^ which 
implies A 2 + 2r]A + fi 2 , = 0, i.e. A = — 77 ± iVL where 



n = yjn 2 - r, 2 . (2.8) 

We assume weak damping, r//0,Q < 1. The two fundamental solutions are to be 
superimposed such that the initial condition (2.7) are satisfied. We obtain 

T(r) = i e _r?r sinftr . (2.9) 

The original equation of motion has thus been solved for an arbitrary oscillating 
electric field strength E = E(t). 

2.3 Susceptibility 

p = — ez is the dipole moment of the electron under consideration. If there are 
N of them per unit volume, the polarization is P = Np. We just have obtained 
the result 

Ne 2 

P(t) = — drr{r)E{t-r) . (2.10) 
™> Jo 

Let us Fourier decompose the electric field strength: 

E{t) = J £ e-™ E{u) . (2.11) 
We work out 

P(t) = eo j ^e~ iut xH^H (2.12) 

where 

, v Ne 2 [°° iur r , , Ne 2 1 

X{lo) = / are T(r) = — . . — . (2.13) 

meo Jo me o + w + ir]){\l — uj — IT)) 

The Fourier component P(uj), which is declared in analogy to (2.11), is propor- 
tional to the Fourier component E{u) at the same angular frequency u: 



P(u) = eo X (co)E(u;) 



(2.14) 
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2 A SIMPLE MODEL 



X = x( w ) is called the susceptibility of the material under study. It is a function 
of angular frequency. Note that we may rewrite (2.13) as 



(2.13) is an oversimplification because not all electrons have the same resonance 
frequency and the same damping constant 77; in fact, one must sum over such 
terms. However, (2.15) or a sum over such terms reflect the essential features 
quite well: 

• The static value %(0) is always positive and real. 

• x( w ) — * with uj — > 00. 

• The imaginary part of the susceptibilty is always non-negative. 

• The imaginary part of the susceptibility is large close to a resonance 



• Dispersion is unavoidable, susceptibilities must be different for different 
angular frequencies. 



xM = 



Ne 2 1 

meo fig — uj 2 — 2irjuj 



(2.15) 
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3 Thermodynamic equilibrium 

In this serious of lectures we will study the interaction of a rapidly oscillating 
electromagnetic field with matter. Therefore, a recollection of basic statistical 
thermodynamics seems to be appropriate. 

3.1 Observables and states 

We describe the system quantum-mechanically, i. e. by an appropriate Hilbert 
space 7i. Observables M, or measurable quantities, are represented by self- 
adjoint linear operators mapping the Hilbert space into itself, M : TC — > TC. 
As is well known, the eigenvectors of self-adjoint operators form a complete 
system of normalized and mutually orthogonal vectors. The eigenspaces of an 
observable M describe the alternatives, the eigenvalues the measured values of 
these alternatives. 

Pure states of the system are represented by wave functions, or vectors of TC. 
For a large system, the notion of a pure state is an oversimplification. In 
fact, there is a complete system of normalized and mutually orthogonal vectors 
01,02, .. . and a set of probabilities w±,W2,-- ■ ■ The Wj are the probabilities 
that the system will be found to be in the pure state (pj. 

There is a unique linear operator W such that the <j>j are its eigenvectors and 
the Wj it eigenvalues, 

Wfy = wrfj . (3.1) 

This probability operator 1 W describes the state of the system. It is character- 
ized by 

< W < 1 and trW = 1 . (3.2) 

Here and 1 stand for the zero and the unity operator, and tr denotes the 
trace. (3.2) boils down to < Wj < 1 and Yl w j = 1- 

The expectation value (M) of the observable M, while the system is in state 
W, is given by 

(M) = ^w j (<f> j ,M<j> j ) = tvWM . (3.3) 

3 

3.2 The first law of thermodynamics 

The energy observable H (the Hamiltonian) is of particular interest because it 
governs the time development of the system. Its expectation value is convention- 
ally denoted by U, the internal energy. It is a convention that thermodynamic 
systems are at rest, hence energy is internal energy. 



1 also density matrix 
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3 THERMODYNAMIC EQUILIBRIUM 



The internal energy of a system may change because the state changes or be- 
cause the Hamiltonian changes, 

SU = tv5WH+ ti W5H . (3.4) 

The first contribution 5Q = trSWH is heat, the second contribution 5 A = 
tr WSH is work 2 . The Hamiltonian usually depends on external parameters 
which may change. An external electric field is a typical example. The state- 
ment SU = 5Q + 5A is the first law of thermodynamics: energy may be trans- 
ferred to a system as heat or as work. 

3.3 Entropy 

A mixed state W is a mixture of pure states. The pure states <pj are contained 
in the mixed state W with probabilities Wj. Two mixed states W\ and W 2 may 
be further mixed to become W = a\Wi + a 2 W 2 where < 01,02 < 1 and 
0:1 + 0:2 = 1. It is a simple exercise to show that < W < 1 holds as well as 
tvW= 1. 

We require a measure to answer the question: how much mixed is a state. This 
measure should vanish if W represents a pure state, and it should increase upon 
mixing. Here we just report the result of a lengthy discussion: 

S(W) = k B J2 w i ln — = -kBtrWlnW (3.5) 

j J 

is a measure of the degree of mixture. S(W) is the entropy of the mixed state 

W. The Boltzmann constant fee shows up for purely historical reasons. 

If all probabilities Wj equal zero, up to one, which must be unity, than the 

entropy vanishes. The entropy of a pure state vanishes. It can also be shown 

that 

aiS(Wi) + a 2 S(W 2 ) < S(aiWi + 02^2) (3.6) 
is true for states W\,W 2 and weights 01,02- Mixing increases entropy. 

3.4 The second law of thermodynamics 

Assume a system which is well isolated from its environment. The internal 
energy within the system remains constant. Nevertheless, the environment will 
influence the time evolution of the system. The second law of thermodynamics 
states that the interactions between the system and its partially chaotic envi- 
ronment will increase the amount of chaos within the system which is measured 
by the entropy of its state. Put otherwise, the entropy of the system's state 
increases in the course of time. 

2 German Arbeit. The symbol W (for work) is already in use (for the probability operator, 
or state) 



3.5 Irreversible processes 
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It increases until it has reached a maximum. The state with maximal entropy 
does not change any more, it describes the equilibrium of the system with 
respect to it environment. Let us denote this equilibrium state by G. Denote 
by 

S = {W:H^H\0< W<l,trW = l} (3.7) 
the set of states. The equilibrium state, or Gibbs state, is characterized by 

S(G) = max S(W) , (3.8) 

which amounts to 

S(G + SW) = S(G) for ti5W = and ti5WH = (3.9) 

for small deviations 5W. G + 5W must be a state, and U = tr (G + 5W)H 

should remain constant. 

The solution of this problem reads 

G=e (F-H)/k B T (310) 

where F and T are Lagrange multipliers. F, the free energy, is a number. H is 
the Hamiltonian of the system, and T the temperature of the Gibbs state. 
The free energy is calculated according to 

F = -k B Tlntr e~ H/kBT (3.11) 
such that tr G = 1 holds. The temperature T is determined by solving 

trHe- H ^ T , , 

for T. One can show that the right hand side of (3.10) increases with the 
temperature T which guarantees a unique solution. Adding energy to a system, 
by heat or work, will make it warmer. 

3.5 Irreversible processes 

We denote the external parameters summarily by A, i. e. H = H(X). The 
free energy depends on the temperature and on the external parameters, F = 
F(T, A). One can easily show that 

dF = -SdT - v rdX r (3.13) 

holds true where 



S = S(G) 



(3.14) 
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3 THERMODYNAMIC EQUILIBRIUM 



is the entropy of the Gibbs state. The generalized forces V r are given by 



We speak of a process if the temperature T = Tt of the system's environment 
and the external parameters A r = A r ,t change in the course of time. The 
process is reversible provided the state Wt of the system is always very close to 
the corresponding equilibrium state, 



If the external parameters change too rapidly, the system will remain away 
from equilibrium, and we speak of irreversible processes. With light, the typical 
spatial dimensions of material points are micrometers, and the corresponding 
period is 3 x 10 14 s. This time is much too short for achieving equilibrium. The 
interaction of light with matter is an irreversible process. We will describe in 
the following sections how to cope with this problem. 



V r = —trG 



dH{\) 



(3.15) 



W t ~G{T u \ t ) 



(3.16) 
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4 Perturbing the equilibrium 

In this section we will describe how the state of a system changes with time. 
We concentrate on small perturbations of the equilibrium state. 



4.1 Time evolution 

So far we have spoken of preparing a state W and immediately measuring an 
observable M. Let us now discuss the situation that we wait for a certain time 
t between preparing W and measuring M. 

Let us first discuss the Heisenberg picture of time evolution. Waiting the time 
span t and then measuring M is a new observable M t . M t must be self-adjoint. 
Its eigenvalues, the possible outcomes of a measurement, are independent on 
the waiting time span t. Therefore, 

M t = U-tMU t (4.1) 



holds where Ut is a unitary operator. 
Because of 



U tl+t2 = U t2 U tl (4.2) 
we conclude that U depends exponentially on t, 

U t = e ~ S tH (4.3) 

where H is a self-adjoint operator, the Hamiltonian, or the energy. Note that 
(4.3) results in 

jM t= '-[H,M t ] , (4.4) 

the Heisenberg equation of motion. [A, B] = AB — BA is the commutator of A 
with B. 

Another aspect is the Schrodinger picture. Conceptually, preparing W and 
waiting a certain time t is preparing a state Wt- Because of 

tr WM t = tr WU-tMUt = tr (,117' h \l (4.5) 

we have to write 

W t = U t WU-t . (4.6) 
We have made use of tr AB = tr BA. 

(4.6) results in the Schrodinger equation for mixed states, namely 

±W,= i[W„H\. (4.7) 
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4 PERTURBING THE EQUILIBRIUM 



4.2 Interaction picture 

Very often the Hamiltonian may be split into a manageable, time independent 
part H and a small perturbation V = Vt which may depend on time. Just think 
of H as the Hamiltonian of matter and Vt describing the perturbation by the 
electric field of a light wave. 

It is then useful to work in the Heisenberg picture with respect to H only. 
In the Schrodinger picture we have to solve 

jW t = l -[W t ,H + V t ] . (4.8) 

With 

A(t) = U-tAUt (4.9) 
where the time translation operator Ut is defined by (4.3) we arrive at 

or 

jW t (t)= l -[ Wt (t),V t (t)} . (4.11) 

This equation of motion is driven by the perturbation only. If there were no 
V = Vt, then Wt{t) would be constant. It is to be expected that a small 
perturbation will cause the state Wt(t) to change only slowly. 

4.3 Perturbing the Gibbs state 

Assume that we perturb the Gibbs state: 

W t -»• G = e ( F ~ H )/ k * T and V t -► for t -»• -oo . (4.12) 

Because the Gibbs state is stationary with respect to H, i. e. G(t) = G, we 
may also write 

W t {t) -> G for t -> -oo . (4.13) 

By combining the initial condition (4.13) with (4.10) we arrive at the following 
integral equation: 

W t {t) = G + f t ds l -[W s {s),V s {s)] . (4.14) 

This integral equation gives rise to a perturbation expansion with respect to V. 
In lowest order we have to set Wt(t) = G + In first order we arrive at 

W t {t) = G+ f ds t[G , V s (s) ] + . . . . (4.15) 

Inserting the first order expression (4.15) into (4.14) will result in the second 
order approximation, and so forth. 



4.4 Time dependent external parameter 
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4.4 Time dependent external parameter 

We now specialize to 

V t = -J2^r(t)V r . (4.16) 

r 

is an external parameter, V r the corresponding generalized force, 
an observable. We may rewrite the linear response (4.15) of the system to a 
perturbation of the equilibrium state into the following form: 

W t {t) = G- f dsYj^s) % -[G, V r (s)} + ... , (4.17) 



or 



W t = G- f dsJ2 X r(s)l[G,V r (s-t)} + ... . (4.18) 

J -co r 

The expectation value of an observable M is 

poo 

trW t M=trGM+ drVA r (t-r)r(MF r ;T) + ... (4.19) 

r 

where 

r(MV r ; t) = tr G l - [M(r), V r ] . (4.20) 

This is a remarkable result. The expectation value of an observable M is its 
equilibrium expectation value plus an additional contribution. The addition 
depends on past perturbations only. It is proportional to the external param- 
eters A r . The Green function r = r(MV r ;r) depends on the age r of the 
perturbation and is linear in the perturbing generalized forces V r . In fact, the 
commutator [M(t),V r (s)] will vanish unless V r (s) and M(t) affect each other. 
What is most important: the Green function is an expectation value of the 
unperturbed state G. Loosely speaking, the Gibbs state G knows already how 
it will react on perturbations. 
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5 DIELECTRIC SUSCEPTIBILITY 



5 Dielectric susceptibility 

We now specialize to ordinary matter in an oscillating electric field. We may 
write 

H = y^ + -^y l %qa , (5.i) 

^ 2m a 4vre \x b - x a \ 

a b>a 

for the Hamiltonian of matter. Charged particles, nuclei and electrons, interact 
via Coulomb forces. a,b . . . label the particles, m a and q a are the respective 
masses and charges. x a and p a denote the particle positions and momenta, 
they are observables. This is no treatise on many particle physics, and (5.1) 
serves only to introduce notation. In particular, spin and relativistic effects are 
missing. 

5.1 Polarization of matter 

If there is an electromagnetic field, the interaction with matter can be classified 
by a multi-pole expansion. For our purpose, the electric dipole approximation 
is sufficient. 
We define by 

P(x) = y j q a x a 5\x-x a ) (5.2) 



the density of dipole moments, or polarization as a field of observables. The 
interaction of the external electromagnetic field with matter is described in 
electric dipole approximation by 

Vt = ~ j d 3 yE(t,y)-P(y) . (5.3) 

We assume a dielectric material, in contrast with a ferroelectric material. With- 
out electric field there is no polarization, tr GP{x) = 0. By inserting (5.3) into 
(4.19) and (4.20) we arrive at the following expression for the time-dependent 
polarization: 

/*oo r 

Pi(t,x) = tr W t Pi(x) = J^ dr J d 3 y r ij (r,x,y)E j (t-T,y) . (5.4) 

Here and later we rely on Einstein's summation convention: a sum over the 
doubly occurring index j is silently understood. 
The tensor of Green functions 3 is given by 

rij{T,x,y) = tTG % -[Pi{T,x),Pj{Q,y)\ . (5.5) 



! Tij(T, x, y) = r(Pi(x)Pj(y); r) in the notation of section 4 



5.2 Dielectric susceptibility 
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Recall the definition of the polarization Pi (x) and its time translation 

Pi(t, x) = U- t Pi{x) U t where U t = e ~ ^ tH . (5.6) 
Also recall the definition of the equilibrium, or Gibbs state 

G=e (F-H)/k B T (5?) 

H in (5.6) and (5.7) is the Hamiltonian (5.1) of unperturbed matter. 

We assume that the Gibbs state is not only invariant with respect to time 

translations, but also with respect to spatial translations. One may then write 

r ij {T,0 = trG l - [^,0,^(0,0)] (5.8) 
instead of (5.5) and 

Pi(t, x) = tr W t Pi(x) = J dr J d 3 d r^(r, £) E 3 (t -t,x~0 (5.9) 

instead of (5.4). 

5.2 Dielectric susceptibility 

Let us Fourier decompose the external electric field: 

Ei(t,x) = J—e J 7^)3 e Ei(uj,q) . (5.10) 

We may likewise decompose the polarisation 4 Pi(t,x). 

The well-know convolution theorem of Fourier theory allows to write 

Pi(uj,q) = e Xij(u,q) Ej(uj,q) (5.11) 

where 



Xij 



(u,q) = - T dr e iUJT [ d 3 £ e iq ^ r tJ (r,0 (5.12) 
eo Jo J 



is the tensor of dielectric susceptibility. It is a simple exercise to verify that the 
susceptibilities are dimensionless. If we could solve all equations which we may 
write down, the susceptibility of any material can be calculated. It depends 
on frequency, on the wave vector, on temperature, and on all other parameters 
which characterize the equilibrium state. 



4 We use the same word for the observable Pi — Pi(x) and its expectation value in the time 
dependent state Wt, i.e. Pi(t,x) = trV^tP^^)- 



20 



5 DIELECTRIC SUSCEPTIBILITY 



5.3 Susceptibility proper and optical activity 

Since the velocity of sound, which sets the time to distance scale in solids, is 
so much smaller than the velocity of light, the susceptibility varies only weakly 
with the wave number. We may expand 



In most optical applications only Xiji^ ) 1S required. The frequency dependent 
tensor Xij( w ) is likewise called the susceptibility of the material in question. 
The next term describes optical activity. Higher order expansion coefficients 
are seldom encountered. Optical activity is deferred to a later section. In 
the following we are concerned entirely with the ordinary frequency dependent 
susceptibility tensor Xiji^)-, 



Xiji^J, q) = Xij{u) + Xijk(v) qk + ■■■ ■ 



(5.13) 




(5.14) 
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6 Dispersion relations 

We will now exploit the fact that the response to a perturbation is retarded. 
The generalized susceptibility is given by 

[•CO 

X (AB;u)= / dr e WT r(AB;r) , (6.1) 
J o 

where the response function is 

r(AB;r) = tvG l -[A(r),B(0)} . (6.2) 

It weighs the impact of a perturbation caused by B on A after a delay r. 
See appendix A for a summary on causal functions. 

6.1 Retarded Green function 

In fact, the generalized susceptibility is the Fourier transform of a product, 

X {AB-uj) = Jdr e luJT e(r)r(AB;r) , (6.3) 

where 6 = 6{x) is the Heaviside 5 jump function, a distribution. 

The Fourier transform of a product is the convolution of the respective Fourier 

transforms, 

X (AB; u) = J ^ 6{u - u) f(AB; u) . (6.4) 
The Fourier transform of the Heaviside function is 

6{u)= r dt e^-^ = . (6.5) 

Jo u + ie v ' 

e is a positive, yet arbitrarily small number. Therefore, 

xiAB-u) = ±-l d ^— f(AB;u) . (6.6) 

2m J u — lo — it 

(6.6) says that the pole at u = uj is to be avoided by running in the lower 
complex plane. 

Now, the Gibbs state is stationary, therefore 

( l -[A(r),B(0)}) = ( l -[A(0),B(-r)}) (6.7) 
holds true which implies 

r(AB;r) = r(BA;-r) (6.8) 



r '6(x) = for x < and 6{x) = 1 for x > 
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6 DISPERSION RELATIONS 



or 



f(AB; uj) = -f(BA; -co) = -f(BA; u>)* . (6.9) 

The last conclusion relies on the fact that response functions r(AB; r) are real. 
With this information we may work out 

X (BA;u,y = ± f ^— r(AB,u) . (6.10) 

2m J u — uj + te 

Note the similarity with (6.6). The only difference is how to evade the singu- 
larity u = lo. The difference of (6.6) and (6.10) is a ring integral around the 
singularity, it results in 

x{AB;lo) - x{BA;lo)* = r(AB;uj) . (6.11) 

We define the refractive part of the susceptibility by 

X '(ABM = X(AB ' M+ 2 X{BA '" Y (6.12) 
and the absorptive part by 

x > lAB;a) = xM»)-«bm u y (6 13) 

Both are Hermitian in the sense that interchanging A and B as well as complex 
conjugating it leave the expressions unchanged. 

Inserting f(AB;w) = 2ix"(AB;uj) into the sum of (6.6) and (6.10) we arrive 
at 

x'(AB;u) = l^J ^ x"(AB;u) . (6.14) 

The principal value of the integral is the mean of avoiding the singularity via 
the upper and the lower complex -u-plane. It may also be defined by 

J U UJ V7-oo Jco+eJ U-UJ 

(6.15) is the prototype of a dispersion relation. In the following subsection we 
will specialize to optics thereby justifying the terminology: refractive, absorp- 
tive, dispersion. 



6.2 Kramers-Kronig relations 
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6.2 Kramers-Kronig relations 

Recall the definition of the dielectric susceptibility tensor 



1 f°° ' f i 

Xij(<") = -J dre~ WT J d^trG-[P t (r,O,Pj(0,0)] ■ (6.16) 

It fits into the general scheme, 

XtfM = j d^x(Pi(OPj(0);r) . (6.17) 

Response functions r(AB; r) and their susceptibilities %(AB; uS) are linear in 
both arguments, A and B, and the same can be said of the dispersion relation 
(6.16). It can be summed over A and remains true. 
With 



and 



X>) = ^^^ 



2i 



we may write 



X» = ipr / -^-X» • (6-20) 

This particular form of a dispersion relation is known as the Kramers-Kronig 
relation. It first of all explains why refraction depends on frequency. A beam 
of natural light consists of many colors, and a glass prism will lead to disper- 
sion. The Kramers-Kronig relation also tells that there is no refraction without 
absorption, although at different frequencies. It remains to show that the ab- 
sorptive part Xiji^) is a non-negative Hermitian tensor. 

6.3 Refraction and absorption 

We discuss Maxwell's equations for purely periodic fields in absence of charges 
and currents: 

V x H = -iuje eE and V x E = i(jfi H . (6.21) 

The remaining divergence equations are automatically fulfilled. 
The permittivity tensor e of (6.21) is 



(6.22) 
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6 DISPERSION RELATIONS 



We simplify the discussion and discuss an isotropic material, €ij(u) = <% e(u>). 
We intend to show that the real part of the susceptibility is responsible for 
refraction while the imaginary part causes absorption. 
Consider a damped plane wave traveling in z direction: 

„ „ — iut ink z —az/2 _ _ 

E\ = Ee e u e ' ; E 2 = ; E 3 = . (6.23) 

The wave number is q = nko where ko = u/c. n is the refractive index, a the 
absorption constant. Inserting (6.23) into (6.21) shows that we have guessed 
correctly, provided 



„ , a \ 2 no 

e = e ' + ie"= [n + i—] = n 2 + % — + .. . (6.24) 
V 2k oJ k 

holds true. Thus we have shown that 

n = y/7' (6.25) 
is the refractive index and 

a = ^ (6.26) 
n 

the absorption constant. Recall that the energy current density S is quadratic 
in the electromagnetic field strengths. Therefore, a is indeed the power decay 
constant of radiation. Proving e" > is close to proving the second law of 
thermodynamics. We will address this challenge in the following section. 
The susceptibility is the Fourier transform of a real function, therefore x( w ) = 
X*(—lo). Consequently, the real part is an even, the imaginary part an odd 
function of angular frequency. Averaging x'( w ) an d x'(~ w ) an d taking into 
account that ux" (u) / (u 2 — to 2 ) is an even function of u we arrive at 

X 'H= 2 Pr f °du U f^ 2 . (6.27) 
6.4 Oscillator strength 

We have derived in section 2 an expression for the susceptibility of elastically 
bound, weakly damped electrons. If there are iV such electrons per unit volume 
resonating at frequency u, the refractive part of the susceptibility is 

„ , Ne 2 1 , 

X '{u) = j g . 6.28 

meo u z — u £ 

This explains why the Kramers-Kronig relation is often written in the following 
form: 

n 2 {io) = 1 + — Pr r du -4^K ■ (6-29) 
me Jo u 2 -uj 2 



6.4 Oscillator strength 
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/ = f(u) is the oscillator strength describing the distribution of resonance 
frequencies. 



N= / duf(u) (6.30) 
Jo 

is the spatial density of elastically bound electrons. The oscillator strength of 
(6.29) is given by 

/H = -^nH«H . (6.31) 
7r e z 

Note that the static dielectric constant is 
Ne 2 1 

e (0) = l + — — (6.32) 
meo uj z 

where to is an average resonance frequency defined by 



i 



30 du J ^= N , . (6.33) 
u z UJ Z 

2 



The integral converges because f(ui) is proportional to ux'{w) oc to 
For very high frequencies (X rays) the refractive index is given by 

n 2 ( W ) = l-— ^ • ( 6 - 34 ) 
meo w 

A refractive index less than 1 means that the phase velocity uj/k = c/n is larger 
than c. This does not imply that X ray signals may travel faster than light. 
Wave packets travel with group velocity v = duj/dk which differs from phase 
velocity. The speed of a signal is yet another story. 

If the angular frequency u lies outside an absorption band, f(u) = 0, the 
principal value operator Pr in (6.29) may be omitted, and we deduce 

dn 2 2uoe 2 f°° , f(u) 

-j— = / du - ' > . 6.35 

du; meo Jo \ u ~ u ) 

The permittivity, or the refractive index, grows with angular frequency. This 
behaviour of dispersion is normal. Abnormal dispersion — the refractive index 
decreases with increasing frequency — is possible only within absorption bands. 
If f(co) > 0, the above argument does not hold true. Abnormal dispersion 
comes necessarily with large absorption. 
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7 DISSIPATION-FLUCTUATION THEOREM 



7 Dissipation-fluctuation theorem 

It is well known that things become simpler if one concentrates on the essentials. 
We shall therefore discuss the relation between dissipation and fluctuation on a 
very abstract level. The interaction of a light wave with matter will then serve 
as an example. 

We discuss a system which, if unperturbed, is described by its Hamiltonian H. 
The system's equilibrium state is G oc exp (—/3H) where (3 = l/k&T is short for 
the inverse temperature (up to the Boltzmann constant &b). The expectation 
value of an observable M in the equilibrium state is denoted by (M) = tr GM 
throughout. 

If the equilibrium is disturbed by a time dependent contribution 



H t = H-X(t)V , (7.1) 
the linear response, as felt by an observable M, is given by 

/•oo 

(M) t = {M)- / dT\(t-T)r(MV;T) . (7.2) 
Jo 

There is an intimate relation between response functions T(AB;t) and corre- 
lation functions K(AB; r) which is the subject of this section. 



7.1 The Wiener-Khinchin theorem 

If A denotes an observable of the system under discussion, 

A(t) = U-tAUt where U t = e ~ ^ tH (7.3) 



is considered to be a process. Note that the equilibrium values are constant, 

(A(t)} = tr GU-tAUt = tr U t GU- t A = ti GA = (A(0)) . (7.4) 

However, there will be fluctuations. We describe them by the time correlation 
function 

K(AB; r) = MC + r)B(t) +B W A(t + t)) _ {A) (J?) ^ 

Because we have chosen the symmetrized product, the correlations function of 
two observables A and B will always be real. Note that the time argument t 
is absent on the left hand side of (7.5). K(AB; 0) is the correlation proper. 
K(AA; t) is a time auto-correlation function. 
Let us insert the Fourier decomposed processes, 

A(t + r) = {A) + J^e-^ t + T U( U ) (7.6) 



7.2 Kubo-Martin-Schwinger formula 
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and 

B(t) = (B) + J ^ BHu') . (7.7) 
There will be a factor e ^ U ^ ^ leading to a dependency on t unless 

In a stationary state, like in the Gibbs state, the Fourier components of fluctu- 
ations with different frequencies are not correlated. Correlated fluctuations of 
different frequencies would cause beating. 

The frequency dependent function S = S(AB; uS) in (7.8) is called a spectral 
density. With it we may write 

K(AB;t)= f — e~ %UJT S{AB-uj) . (7.9) 

J 27T 

5(^4^4; oj), the spectral density of an auto-correlation function K(AA; r) is never 
negative, 

K(AA;t) = [ — e~ lUJT S(AA;lo) where S(AA;to) > . (7.10) 
J 2ir 

This finding is known as the Wiener-Khintchin theorem. 

7.2 Kubo-Martin-Schwinger formula 

Let us define 

A{z)=e h ztl Ae h zn (7.H) 
for zGC. 

We now exploit the fact that both the Gibbs state and the time translation 
operators are exponential functions of the energy: 

A(z)e- (3H = e~ PH e (3H A(z)e- (3H , (7.12) 

i.e. A(z)G = GA(z — ih/3). We multiply from the right by B and apply the 
trace operator. The result 

(BA(z)} = (A{z - ih0)B) (7.13) 

is the famous Kubo-Martin-Schwinger, or KMS formula. 
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7 DISSIPATION-FLUCTUATION THEOREM 



7.3 Response and correlation 

Note the formal similarity between the response function 

r(AB;r) = ( l -[A(r),B}) (7.14) 
and the correlation function 

K(AB;T)=W T > B t BAiT)) -{A){B) . (7.15) 



The difference is essentially the difference between a commutator and an anti- 
commutator. 

We therefore investigate the product, 

4>{t) = {A{t)B)-{A){B) . (7.16) 
Its Fourier transform is 



(u) = / dr e lUJT <f>(r) . (7.17) 



Let us now discuss the function 

f(z) = f < £e- iuz fa) . (7.18) 

It can be shown that this function is analytic within a large enough strip of 
the complex z-plane around the real axis. On the real axis we clearly have 
/(r) = 4>(t). Indeed, function / is the analytic continuation into the complex 
plane of function (j). 
With (7.11) we may work out 

g(z) = {A(z)B)-{A)(B) (7.19) 

which is analytic in a sufficiently broad strip around the real axis. On the real 
axis we have g{r) = /(t). Therefore, g{z) and f(z) coincide everywhere. 
Applying the KMS formula yields 

{BA(r))-(A){B) = f(r-ih(3) . (7.20) 

One should compare this with 



(A(r)B)-(A)(B) = f(r) 



(7.21) 
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7.4 The Callen-Welton theorem 

The response function can now be written as 

r(AB;r)= {/(r) - f(r - ih/3)} . (7.22) 
For the time correlation function we obtain 

K(AB;r) = ^{f(r) + f(T-ih/3)} . (7.23) 
We insert the representation (7.18) of /, 



(7.24) 



K(AB; r) = \\%e~^ 0(c) { 1 + e } . (7.25) 

One can read off immediately the Fourier transforms, and eliminating <fi results 
in 

r(AB;cj) = 2 ^S(AB;u)tanh^- . (7.26) 

lb _■ 

We are more interested in the (generalized) susceptibility 

j(AB;u)= dre lUJT r(AB-T) . (7.27) 
Jo 

Recall the result 

7(AB;w) -j(BA;u)* = r(AB;u) (7.28) 

of section 6 where A was Pi and B = Pj. 
Combining (7.26) and (7.28) results in 

7(AB=")-7(B*")- =S( ^. M) l^g^ . (729) 

This is the famous fluctuation-dissipation theorem of Callen and Welton. 
7.5 Energy dissipation 

Let us discuss a process Wt with the initial condition W-oo = G oc e 
This process is driven by the time dependent Hamiltonian 

H t = H-Y J K{t)V r , (7.30) 
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and we assume X r (t) — ► for \t\ — > oo. During the small time span eft the 
following amount of work 

dA = - ^2d\ r (V r ) t = - ^2dt\ r (t)(V r ) t (7.31) 

is spent on the system, where (V r )t = tr WtV r . This follows from the first law 
of thermodynamics. The total work spent on the system is 

A = -J2 f dtX r (t)(V r ) t . (7.32) 

r 

We assume that the external parameters are always small, so that the linear 
approximation to the response is sufficient: 

poo 

(V r ) t = {V r ) + y2 dr r rs (r) X s (t - t) . (7.33) 
s J ° 

The response functions are 

r rs (r) = r{V r V s -T) = ( l -[V r (r), V s (0)}} . (7.34) 
The constant terms (V r ) do not contribute, and we may write 

A = - [ dt I dr X r (t)X s (t - r) 0(t) T rs (r) . (7.35) 

rs J 

This is an expression of type 

J dt' J dt" f(t')g(t'-t")h(t") = J ^j(u)*g(u)h(u) , (7.36) 
and we may write 

A = -^2 j ^ iu X r (u)* Xrs(v) X s (uj) . (7.37) 

rs 

Note that 

,(w) = J dr e iUJT 0(t).T„(t) (7.38) 



X 



are generalized susceptibilities. 

Complex conjugating (7.37) and interchanging the summation indices r, s will 
not change the value of A, 



A 

sr 



£ J ^,A s Hx;HA;H , (7.39) 



7.5 Energy dissipation 
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and averaging (7.37) and (7.39) results in 




ooduj <- 



\ r {uj) 



2i 



(7.40) 



With the dissipation-fluctuation theorem (7.29) we finally arrive at 



A = 




q{(3hw) ^2 A r (u;)*S' rs (a;)A s (u;) 



(7.41) 



where 



q{x) = — tanh 



(7.42) 



is a quantum mechanical correction. By the Wiener-Khintchin theorem the 
matrix S rs = S ts {uj) of spectral densities is non-negative. 

We conclude that it is impossible to perturb an equilibrium state in such a way 
that the system extracts heat from its environment and delivers work: A > 0. 
This statement is an alternative formulation of the second law of thermody- 
namics. We have proven it here for small perturbations of the equilibrium by 
rapidly oscillating external parameters. 

A more prosaic formulation reads: the absorptive part of the susceptibility truly 
causes absorption of energy. 
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8 ONSAGER RELATIONS 



8 Onsager relations 

If the equilibrium is perturbed by more than one driving force, the generalized 
susceptibilities Xrs(^) form a matrix. We have shown in section 7 that the 
imaginary, or absorptive part Xrs{^) — Xsr{^)* is non-negative and proportional 
to the matrix of spectral densities. This was the subject of the dissipation- 
fluctuation theorem. 

In this section we discuss symmetry properties of the susceptibility matrix. 



8.1 Symmetry of static susceptibilities 

Let us first discuss the symmetry properties for adiabatic processes. Let A = 
Ai, A2, . . . denote external parameters or deviations thereof such that 



H(X) = H-J2\ r V r (8.1) 

r 

is the relevant Hamiltonian. If the external parameters change slowly, the state 
W of the system will always be the equilibrium state to the current external 
parameters A. The free energy is given by 

F(X, /?) = - * In tr e ~^ H ^ . (8.2) 

With (...) as expectation value in the Gibbs state to H(0) and (. . .)\ to H(X) 
we find 



(Vr)x=(V r ) + ^Xs d{Vr)X 



dX 



(8.3) 

A=0 



Comparing with (7.33) we conclude 

x „ ( o)- J f*r f .(r)--$gM . (8.4) 

It follows immediately that 

Xrs(O) = Xsr(O) (8.5) 

holds true. Can this symmetry be generalized to fast processes? Almost, as we 
shall see. 

8.2 Time reversal 

To each linear operator A we assign a time-reversed operator A* by 



{Acf))* = 



(8.6) 



8.2 Time reversal 
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where <f> = 4>{x\,X2, ■ ■ • ) is an arbitrary wave function depending on the positions 
of the particles. Let us ignore spin at this point. As usual, * denotes complex 
conjugation. 

Note that the mapping A — > A* is anti-linear in the sense that A + B map to 
A* + B*, but zA becomes z*A* 

A position operator x amounts to the multiplication of the wave function by 
its argument, which is real. Therefore, x* = x. Momentum operators arc 
represented by p = —ih'V, therefore p* = —p. This is the reason why the time- 
reversal is sometimes called momentum reversal. By the way, angular momenta 
change their sign at timer reversal as well. 

If (j> is an eigenfunction of A with eigenvalue a, then (jf is an eigenfunction of A* 
with eigenvalue a*. Hence, if M is an observable, so is M*, and if W represents 
a mixed state, so does W*. Observables are characterized by real eigenvalues. 
States also have real eigenvalues which are probabilities. It is a simple exercise 
to prove the expectation value of an observable in a state coincides with the 
expectation value of the time reversed observable in the time reversed state, 

ti WM = tr W*M* . (8.7) 

An observable is of parity rj (with respect to time reversal) if M* = rjM holds 
true. Because of M** = M only i] = ±1 are possible. Position operators have 
even, momentum or angular momentum operators odd parity. 
Ordinarily, the Hamiltonian has even parity. However, if there is an external 
quasi-static magnetic induction field B, there is a contribution — ^ fJ-iBi, the 
magnetic dipole interaction. Since the magnetic momentum is proportional to 
an angular momentum, which has odd time reversal parity, we conclude 

H(B)* = H{-B) . (8.8) 

The external magnetic induction field has to be reversed in order to guarantee 
even parity of the Hamiltonian. 

Now, the Gibbs state depends, among other external parameters, upon the 
quasi-static magnetic induction applied to the system, G = G(B). We conclude 

G{B)* = G{-B) . (8.9) 

The time translation operator is 

U t (B) = e~X tH W , (8.10) 

and we easily work out 

U t (B)* = U-t(-B) . (8.11) 

This explains why we speak of time reversal. The time reversed time translation 
operator U£ translates by the negative time span, U-t- However, B — ► — B has 
to be observed. 
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8.3 Onsager theorem 

Let us come back to the original situation that there is a Hamiltonian with a 
time-dependent perturbation, 

H t = H(B)-J2*s(t)V s . (8.12) 

s 

For the expectation values of the perturbing observables we have found 

/•oo 

trW t V r = tvG(B)V r + / dr r rs (B;T)X s (t-T) , (8.13) 

Jo 

where 

r rs {B- T ) = ti G(B)^[U- T (B)V r U T (B),V s ] . (8.14) 
The latter expression may be rewritten as 

tr G{Bf % - [ U- T (B)V r U T (B) , V s }* (8.15) 
which becomes 

-tvG(—B) ^ [U T (—B)V*U- T (B) , V*] . (8.16) 

We assume that the V r observables have parities rj r with respect to time reversal. 
Thus we may write 

-rj r r] s trG(-B) % - [ U T (-B)V r U- T (B) , V s } (8.17) 

holds. By shifting the entire commutator in time by r (which is allowed because 
G is stationary) we arrive at 

- VrVs tr G(-B) % - [ V r (B) , U^ T V S U T (-B) } , (8.18) 

and reverting the commutator finally yields 

r rs (B;T)=r ]r r ]s r sr (-B;r) . (8.19) 

This symmetry is maintained if (8.19) is Fourier transformed over positive time 
spans r resulting in a symmetry relation for the (generalized) susceptibilities: 

Xrs(B;^) =VrVsXsr{-B;uj) . (8.20) 

(8.20) is an Onsager relation. It applies in particular to the dielectric suscepti- 
bility where rj r = r] s = 1. 

By the way, (8.20) and (8.5) do not contradict each other. r/ r / r] s will imply 
Xrs(0;0) = 0. 



8.4 Onsager relation for kinetic coefficients 
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8.4 Onsager relation for kinetic coefficients 

Let us define the flux of V by 

J = V= l -[H, V] . (8.21) 
In the case of electric polarization P r (x), 

J r = J2laX a S 3 (x-x a ) (8.22) 

a 

are the three components of the electric current density. 

If V has parity r] with respect to time reversal, then J has parity —rj. 

Because the Gibbs state is stationary, the equilibrium expectation value of fluxes 

vanishes. The expectation value J(t) = trWtJ in the perturbed state is also 

called a flux. Fluxes exist because the equilibrium is perturbed. 

With 

r{J r V s -T) = ( \[Vr{r) ,V S }) = r sr (r) (8.23) 
we conclude 

POD 

Mt) = Y, dT r rs (T)\ s (t - t) . (8.24) 

s J 

Fourier transforming results in 

J r (u) = '^2a rs (uj)X s (uj) , (8.25) 

s 

where 

poo 

a rs {uj)= dr e luJT r rs (r) (8.26) 
Jo 

are kinetic coefficients. Differentiating (8.19) and inserting into (8.26) leads to 

a rs (B]io) = r] r r] s a sr (-B;uj) . (8.27) 

This symmetry relation for kinetic coefficients was derived already in 1931 by 
Onsager, although for stationary currents only (a; = 0). His reasoning was 
different and, from today's point of view, questionable. However, Onsager also 
considered other causes of irreversible effects like gradients of temperature and 
chemical potential. 
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8.5 Electrical conductivity and Hall effect 

If the equilibrium state of matter is disturbed by an electric field, then 

j i (uj) = ^G ij (B-u)E j (u) . (8.28) 

3 

This is Ohm's law. <jjj is the conductivity tensor which may depend on an ex- 
ternal magnetic field. For a sufficiently weak magnetic field and for an isotropic 
material we may expand according to 

aij(B; to) = a(u) <% + h(u) ^ e ijk B k + ... . (8.29) 

k 

The first term describes electric conduction in the absence of a magnetic field. 
The second term must be antisymmetric in i and j because it is linear in 
B. There is a current contribution perpendicular to the driving electric field 
strength E and perpendicular to the quasi-static induction B. 
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9 Electro- and magnetooptic effects 

Let us recall the expression 

Xij(u) = -J dre luJT J d^tvG- [ P,(r, £) , P,(0, 0) ] (9.1) 
for the dielectric susceptibility tensor. 

Pi(r,0 = ^qaXaS 3 (xa-Z) (9.2) 

a 

is the polarization at location translated by a time span r. The suscepti- 
bilities are integrals over expectation values in the Gibbs, or equilibrium state. 
Therefore, they depend on all parameters which describe the thermal equilib- 
rium state: composition of matter, temperature, external quasi-static electric 
and magnetic fields, stress, and so forth. 

We will concentrate on quasi-static external electric and magnetic fields E and 
B, respectively. 

A systematic power series expansion up to second order reads 

Xij(E,B) = X °° (9.3) 

+ x\%E k (9.4) 

+ X-} k B k (9.5) 

+ XwEkEt (9.6) 

+ xljkeEkBi (9.7) 

+ XwB k B t (9.8) 

+ ... 

Recall that we speak of the Hermitian refractive part, 

Xij (E,B) = Xii( E > B ) > ( 9 - 9 ) 
and that Onsager's relation demands 

Xij{E,B) = X ji{E,-B) . (9.10) 

The same is true for the permittivity tensor = Sij + Xij ■ 

9.1 Crystal optics 

We assume a real symmetric, but otherwise arbitrary permittivity tensor. The 
medium is assumed to be non-magnetic. Let us look for plane wave solutions: 

F(t,x) = f e e u . (9.11) 
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9 ELECTRO- AND MAGNETOOPTIC EFFECTS 



F is any vector field, w, a unit vector, defines the direction of propagation, ko 
is short for uj/c, and n denotes the refractive index. 
Maxwell's equations read 



cb 



nwx(ee) ; e = nwx(cb) 



(9.12) 



Note that the remaining equations w ■ (e e) = as well as w ■ b = are 
automatically fulfilled: the dielectric displacement D and the induction field B 
are divergence free, c is the speed of light in vacuum. 

For a prescribed direction w of propagation, (8.2) may be understood as an 
eigenvalue problem, w x ... is a linear operator described by the matrix 
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(9.13) 



Eliminating the induction result in 



e~ 1 W 2 e 



n 2 e 



(9.14) 



Because of W w = 0, there is one unphysical solution, namely e = w, with 
n~ 2 = 0. The remaining two eigenvectors describe the polarization of allowed 
propagation modes, the eigenvalues rC 2 defining the respective refractive in- 
dices. 



Optically isotropic media The three eigenvalues of the susceptibility coin- 
cide, Xij = X$ij- Amorphous substances, such as glass, and crystals with cubic 
symmetry, such as NaCl, are examples. All transversely polarized plain waves 
are eigenmodes, the refractive index being n = y/1 + \. 

Optically uniaxial crystals If two eigenvalues of e^P coincide, but differ 
from the third, one has a preferred axis and, orthogonal to it, a preferred plane. 
If the polarization vector lies in the preferred plane, the beam propagates with 
the ordinary main refractive index 6 n°, the square root of the doubly occurring 
eigenvalue. If the beam is polarized in the preferred axis, one speaks of the 
extraordinary main refractive index n e . The preferred axis is an optical axis in 
the following sense: if light propagates along it, the refractive index does not 
depend on polarization. 

If a light beam enters an optically uniaxial crystal at an arbitrary angle, it 
splits into an ordinary 7 and an extraordinary beam which travel with different 
refractive indices. One therefore speaks of double refraction, or birefringence. 
LiNb03 is a well-known and extensively studied birefringent crystal. 

6 The refractive index n is a property of the wave. With fco = ui/c and w as propagation 
direction, the wave vector is q — nkow. 

7 ordinary, because it propagates in the plane of incidence spanned by the incoming beam 
and the surface normal 



9.2 Pockels effect 



39 



Optically biaxial crystals If the symmetry of a crystal is low enough, all 
three eigenvalues of the susceptibilty tensor Xij may be different. One can, 
however, demonstrate that there are two optical axes such that beams, when 
traveling along it, have a polarization independent refractive index. Therefore, 
crystals with three different susceptibilty eigenvalues are called optically biaxial. 
They are birefringent as well. An incident beam splits into two beams which 
propagate with different refractive indices. However, in the general case, none 
of them is ordinary in the sense that it remains in the plane of incidence. 

9.2 Pockels effect 

The addition xjjk^k is linear in the strength of an external electric field. (9.9) 
and (9.10) require xjjk t° be a real third-rank tensor which has to be symmetric 
in the first two indices. 

Note that space inversion t — > t and x — > —x is a symmetry of Maxwell's 
equations if E — > — E and B — ► B. Put otherwise, the electric field strength 
is a polar, the magnetic field strength an axial vector. If the Gibbs state is 
invariant with respect to space inversion, xjjk wm acquire a minus sign, it must 
coincide with its negative, hence vanish. We conclude that a third rank tensor 
is possible only for crystals which do not have an inversion center. If there is 
a location (at x = 0) such that x — ► —x is a symmetry, then xjjk necessarily 
vanishes. 

Lithium niobate 8 LiNb03 may serve as an example. There is an xy-plane with 
a 120° symmetry and mirror symmetry x — > — x. The c-axis, orthogonal to 
this plane, has a preferred direction. After all, LiNb03 is ferro-electric. Conse- 
quently, there are four different third-rank tensors which fulfill all requirements. 
We have deferred a detailed discussion to appendix B. 

The Pockels effect, the dependence of the refractive index on the first power of 
a quasi-static electric field strength, allows to switch and modulate light. 

9.3 Faraday effect 

Faraday has discovered that a magnetic field may affect the propagation of 
light. The contribution X^jk^k to the susceptibility must be Hermitian and 
symmetric if B is reversed. Consequently, Xijk 1S to be purely imaginary and 
antisymmetric in i and j. 

To be specific, let us discuss yttrium iron garnet (YIG), a complicated artificial 
crystal which is transparent in the micrometer wavelength region and ferri- 
magnetic at the same time. There is a contribution 

X °} k B k = iKe i3k M k (9.15) 

to the susceptibility. Because YIG is ferri-magnetic, it is customary to refer to 
the magnetization M instead of B. 



at room temperature of 3m-symmetry 
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9 ELECTRO- AND MAGNETOOPTIC EFFECTS 



Let us discuss a situation where a wave propagates along the direction of mag- 
netization, the z-axis, say. With ko = uj/c and k = n/co(0,0,l) we have to 
solve 



-k x k x E = kieE 



where 



e = 



V 



e iKM 
-iKM e 
e 



/ 



(9.16) 



(9.17) 



With 



eL 



V2 



v / 



and eR 



V2 



( 1 \ 

— i 

V J 



(9.18) 



we have made out solutions which describe left and right hand circularly polar- 
ized light, 



-C'L.R — -fce LjR e e 
Their refractive indices differ: 

nl = e + ifM and n| = e - KM 



(9.19) 



(9.20) 



If a linearly polarized plane wave enters a crystal with linear magnetooptic 
effect, it splits into left and right hand circularly polarized waves. These waves 
travel with slightly different phase velocity. When exiting the crystal, the two 
waves recombine again to a linearly polarized wave. However, the direction of 
polarization is rotated by an angle Q£, where t is the path length through the 
magnetooptic medium. The specific Faraday rotation constant is given by 



e = 



2-7T KM 
T In 



(9.21) 



where n = \fe is the average refractive index. 

The Faraday effect distinguishes between forward and backward propagation of 
light and allows to build an optical isolator. 



9.4 Kerr effect 

Unlike the Pockels effect, which is linear in the external electric field strength 
and requires a medium without inversion center, the Kerr effect is usually much 
weaker, because of second order. It is a property even of isotropic media since 

Xgk oc SijSu (9.22) 

fulfills all symmetry requirements: real, symmetric in the first pair of indices, 
symmetric in the second pair of indices. (9.18) amounts to a shift of the refrac- 
tive index by a term proportional to \E\ 2 . 



9.5 Magneto-electric effect 
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9.5 Magneto-electric effect 

There is no accepted name for a contribution xjjke^kBe to the susceptibility- 
tensor. One may speak of a linear dependence of the Pockels effect on an 
external magnetic field. One may likewise say that the Faraday effect depends 
linearly on an external electric field. The x\}kt mus t be purely imaginary and 
antisymmetric in the first pair of indices. The crystal should be ferro- or ferri- 
magnetic, transparent, and must not have an inversion center. 

9.6 Cotton-Mouton effect 

The Cotton-Mouton effect, a permittivity shift which is of second order in 
the magnetic field, is rather weak for para-magnetic material. However, if it 
occurs in a ferri- or ferro-magnetic substance, it may be strong. A contribution 
Xijke ^ fiij^ke wu l change e by a term which is proportional to \M\ 2 and will 
escape attention. The Cotton-Mouton effect is reciprocal, unlike the Faraday 
effect it does not distinguish between forward and backward propagation. 
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10 SPATIAL DISPERSION 



10 Spatial dispersion 

Recall the definition of the susceptibility tensor: 

Xij (u,q) = - T dr e iuJT [ d 3 £ e ~ iq ^ A;(r,£) (10.1) 

where 

r l3 { T ,i) = tvG 1 - [Pi(T,i), Pj(o,o)} . (io.2) 

The polarization — as an operator — is defined by 

Pi(x) = Y,QaXi5 3 (x a -x) . (10.3) 

a 

The sum extends over all charged particles, and x a is particle's a position, 
likewise an observable. The time shift is accomplished by the unperturbed 
Hamiltonian H. The entire Hamiltonian is 

H t = H- j d 3 xP j (x)E j (t, x) (10.4) 

where the external electric field may depend on time. 
The response functions Pij appear in 

tr W t Pi(x) = J drj d 3 C ^{t, Ej{t -t,x-£) (10.5) 

where W = Wt is the disturbed state, and the right hand side refers to the 
lowest order approximation 9 . 

Fourier transforming E = E(t, x) and the left hand side of (10.5) results in 

P%{w,q) = Xij(u,q)Ej(uj,q) . (10.6) 

10.1 Dispersion relation 

So far we have argued that the dependence of Xij^iQ) on the wave vector 
q is so weak that the value at q = suffices. Interactions in a solid spread 
out with the speed of sound while the electromagnetic field is governed by the 
speed of light, and both differ by orders of magnitude. Still, there are certain 
effects where spatial dispersion — the dependency of the susceptibility on the 
wave vector — has to be taken into account. 

The arguments u) and q for E = E(u,q) are not independent. After all, the 
electric field has to obey Maxwell's equations which, for this purpose, read 

(q 2 5ij - qiq^Ej = {^) 2 (Sij + Xij)Ej • (10.7) 



9 Ferro-electric materials contain an additional static contribution P.* = tr GPi 



10.2 Optical activity 
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Recall that both Ej and Xij m (10-7) depend on w and q. The wave equation 
may be viewed as an eigenvalue equation: q given, find an to such that (10.7) 
is satisfied for a non-trivial eigenvector E. 

The relation between wave vectors q and angular frequency uj is a dispersion 
relation. Since w enters (10.7) twice, as a factor on the right hand side and as 
an argument of the susceptibility, the eigenvalue problem is non-linear. 
This complication is seldom realized because the approximation Xijif^iQ) ~ 
Xiji^i 0) is used. The angular frequency then may be considered as an indepen- 
dent variable, and the refractive index squared 10 appears to be the eigenvalue 
to be determined. 

10.2 Optical activity 

To take spatial dispersion into account, let us expand the susceptibility as 

Xij (uj,q) = Xi j (u) + X™ k qk + -.. • (10.8) 

The additional contribution A%°? = XijkQk causes optical activity by traditional 
terminology. 

It is not difficult to show that invariance with respect to time reversal 11 results 
in 

r ij (T,$) = r ji (T,-t) , (io.9) 

which, in turn, implies 

Xij(uj, q) = Xji(u, -q) ■ (10.10) 

We conclude that Xijk 1S purely imaginary and antisymmetric in the first two 
indices, because the susceptibilty must be Hermitian and because of (10.10). 
We write 12 

&Xij = itijk 9k with g k = G ke qe . (10.11) 

g is the gyration vector which depends linearly on the wave vector. 
With respect to space inversion, G k e must be a rank 2 pseudo-tensor because 
€ijk is a rank 3 pseudo-tensor while the susceptibility and the wave vector are 
normal tensors of rank 2 and 1 , respectively. 

Only crystals with a built-in screw sense can be optically active. Quartz is an 
example, or a solution of natural dextrose. Note that the expression 

Axfj = iKe ijk M k , (10.12) 



10 Recall \q\ = nuj/c where n is the refractive index. 
11 in the absence of magnetic fields 

12 eijk is the total antisymmetric Levi-Civita symbol, a rank 3 pseudo-tensor with respect 
to space inversion 
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10 SPATIAL DISPERSION 



which causes the Faraday effect, has the same structure. In full generality, 
also the Faraday effect should be described by a gyration vector which depends 
linearly on the magnetization vector. However, since the magnetization is a 
pseudo-vector, the tensor G in = G^Mn is an ordinary rank 2 tensor, such 
as fe- 
Optical activity results in the rotation of the polarization vector, either to the 
left or to the right 13 . Also K in (10.12) can be positive or negative. The princi- 
pal difference between the Faraday effect and optical activity is the following. If 
a wave passes the medium twice, forward and backward, the Faraday rotations 
add, the rotation of the polarization vector, as caused by optical activity, is 
reverted. Optical activity is a reciprocal effect. 

The fact that all plants produce glucose in its optically right active form (dex- 
trose) is very astonishing. Artificial glucose consists of equal amounts of left and 
right optically active molecules. One molecule is the mirror image of the other 
form, and no chemical mechanism is known which prefers right-handedness over 
left-handedness. May it be that all plants stem from one and the same mother 
plant? 

As for quartz, both forms occur in nature. The crystal growth mechanism is 
such that entire crystals are optically either left or right active (although there 
are twins). 



Right is clockwise as seen by an observer facing the light source. 
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11 Non-linear response 

Denote by H the Hamiltonian of the unperturbed system. The time translation 
operator associated with it is 

U t = e~^ tH , (11.1) 
as derived in section 4. A time translated operator A is denoted by 

A(t) = U-tAUt ■ (11.2) 
We want to solve the following Schrodinger equation for the system's state Wt- 

jW t = l -[Wt,H + V t ] . (11.3) 

Vt is a possibly explicitly time dependent perturbation of the system. Before 
the perturbations has been switched on, the system was in an equilibrium state, 

W t -»• G = e P( F ~ H ^ for t-» -oo . (11.4) 

Recall that F is the system's free energy which guarantees trG = 1, and 
(3 = where T is the temperature of the Gibbs state. 

As explained in detail in section 4, (11.3) is best solved in the interaction picture, 
with states and observables time translated according to (11.2). One arrives at 
the following integral equation which combines the differential equation (11.3) 
and the initial condition (11.4): 

W t (t) = G + ^ ds l -[W s (s) , V s (s)} . (11.5) 

A sensible result: the state now, at time t, depends on previous states and on 
previous perturbations only, and on the initial state. Small perturbations will 
cause only small deviations from the initial state. 

11.1 Higher order response 

The zeroth order approximation is Wt(t) = G+ The first order contribution 

is obtained by inserting the lowest order into (11.5): 



/_ 



oo 



ds l -[G,V s (s)\ + ... . (11.6) 



The second order addition is 

rt rs 



J ds j S du l -[ l -[G, V u (u)}, V s (s)} + ... . (11.7) 
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The expectation value of an observable M is trGM(t) = tiGM in lowest 
approximation. The next term (linear response) is 

/ dstvG-[V s (s), M(t)} . (11.8) 

J -co h 

The quadratic response may be written as 

f ds J" dutvG l -[V u (u), l -[V s (s), M(t)}} . (11.9) 

Let us now specialize to 

V t = -X(t)V . (11.10) 
The linear response contribution (11.8) is 

poo 

/ dnXit-r^r^in) (li.n) 
Jo 

with 

r( 1 )(r 1 )= tTG^[M,V(-n)] . (11.12) 
(11.9) reads 

POO /*00 

/ dTl / dT 2 X(t-T 1 )X(t-T 1 -T 2 )r^(T 1 ,T 2 ) (11.13) 

with 

r (2) (ri,r 2 )= trG^[^[M,F(-ri)],y(-Ti-T 2 )] . (11.14) 

Note that also the second order response function r^ 2 \ri, t 2 ) is an expectation 
value in the unperturbed Gibbs state. 

11.2 Susceptibilities 

Let us Fourier transform m(t) = tr WfM — trGM and X(t). We arrive at the 
following expression, 

mH = X (1) MAM + J ^- x ( 2 \ u , UJ -u)X(u)X(u J -u) + ... , (11.15) 



with 



/>oo 

x«M = / dne^VWln) , (11.16) 

JO 



1 1.3 Second harmonic generation 
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X (2) (-i,-2)= rdne^ + ^i r dr 2 e W ^ r^{r^r 2 ) ,(11.17) 
Jo Jo 

and so forth. 

For the interaction of the electromagnetic field with matter, in electric dipole 
approximation, we replace (11.10) by 

V t = - J d 3 xE,(t,x)P t (x) . (11.18) 

E(t, x) is a time dependent external electric field, P(x) the polarization at x. 
As usual, we use the same symbol for the polarization (as an observable) and 
its expectation value. It Fourier transform is 

P i {u) = X ^\u)E j {u) + J ^ X $(u,u-u)E j (u)E k (u-u) + ... (11.19) 

The expression 

xgVi) = jT J dr ie lUJlTl Jd^( l - [P,(0,0) , ^-(-Ti.O ]> (H-20) 

for the linear response susceptibility has been thoroughly studied in previous 
sections. 

The quadratic response is described by 



XiJfc(wi,W2) = 

J™ d Tl e 1 ^ + J™ dr 2 e^ T2 J d 3 £ J d 3 r, ("-21) 

( ^ [ l - [P,(0,0) , ^(-T!,^)] , i\(-ri - r 2 ,77)]) . 



(11.21) is already a rather complicated expression. What is the counterpart 
to the Kramers-Kronig relation? What is the counterpart to the fluctuation- 
dissipation theorem? Does the second law of thermodynamics hold up to second 
order response theory? Are there symmetries with respect to frequencies and 
tensorial indices? Non-linear response theory is not yet well investigated, and 
we shall stop at this point. 



11.3 Second harmonic generation 

(2) 

X i j k { U} ii ijJ 2) is a proper tensor of rank three, symmetric in the second and 
third index. It vanishes if the Gibbs state is invariant with respect to space 
inversion. Only crystals without inversion center will respond in second order 
to perturbations by an electric field 14 . 

14 Xijfc(0, u>) is the same as xfjki 10 ) which describes the Pockels effect. 
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11 NON-LINEAR RESPONSE 



Let us discuss lithium niobate, LiNbC>3. We choose a coordinate system such 
that the z-axis coincides with the crystallographic c-axis. See appendix B for 
details. Assume a light wave traveling with wave number q in x-direction, being 
polarized in y-direction. The relation between angular frequency u> and wave 
number is 

q = n°(u;)" , (11.22) 
where n° is the ordinary refractive index. The plane wave is described by 

E 1 = E 3 = and E 2 (t,x) = A e~ lUjt e lQX . (11.23) 

(2) (2) 

X122 an d X222 van i sn ) therefore the second order polarization response is 

P 1 = P 2 = and P 3 (t, x) = A 2 xg } 2 e ~ 2lUJ e 2iqX . (11.24) 

In general, this does not excite a light wave because frequency and wave number 
do not fit. Only if the phase matching condition 

2q = n c (2uj) 2 ^- (11.25) 

holds, (11.24) will excite an electromagnetic plane wave. Note that it propagates 
with the extraordinary refractive index. 
(11.22) and (11.25) amount to 

n°(u) = n c {2u) . (11.26) 

The ordinary refractive index of lithium niobate is always larger then the ex- 
traordinary. Since both increase — in the infrared or visible — with frequency, 
condition (11.26) can be met, for a certain angular frequency ui. This frequency 
to be doubled depends on all quantities which affect the refractive indices, such 
as temperature, composition of the material etc. 

Note that the intensity of frequency doubled light grows quadratically with the 
intensity of incident light. In particular, a resonator which confines w-light and 
is leaky for 2w-light, allows high conversion rates. 

Frequency doubling or tripling is technically important because cheap semi- 
conductor lasers emit light of low frequency, while many optical applications 
demand short wavelengths. 

Frequency doubling is well known in music. Musical instruments are highly non- 
linear. A harmonic tone (of one frequency only) always excites other tones, such 
as the octave, or the second harmonic. Therefore, frequency doubling is also 
known as second harmonic generation, or SHG. 
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A Causal functions 

Assume 

/w = /^" !W Sh (a.i) 

where g = g(u) is analytic. It is a simple exercise to prove 

f(t) =0 for t < (A.2) 

provided g = g(u>) is holomorphic 15 in the upper plane Im u> > 0. The converse 
is also true. If / = f(t) is a causal function, as characterized by (A.2), its 
Fourier transform g = g(co) is holomorphic in the upper half plane. 
Hence, f(t) = 9(t)f(t) holds true which may be expressed as 

J 2it e — i{uj — u) 

where the limit < e — ► is understood. The Fourier transform of a product 
is the convolution of the respective Fourier transforms. 

(A. 3) says that the pole at u = to has to be avoided in the lower plane. Adding 
and subtracting an integration path which avoids the singularity u = to in the 
upper half plane results in 

/x ^ I du q(u) , . 

g(u ) = -iPi / — . (A.4) 

Pr J denotes the principal value integral, the average of avoiding the singularity 
in the lower and the upper half plane. The dispersion relation 

// \ ^ f du g"(u) . . 

g'(to) = Pr / — ^L^J. (A.5) 

is a simple consequence, where g 1 is the real and g" the imaginary part of g. The 
Fourier transform of a causal function obeys a Kramer-Kronig like dispersion 
relation. 

So far we never required / = f(t) to be real. If / is a causal function, then if 
is causal as well. Therefore, 

is an alternative formulation. 

If / is also real, then its Fourier transform must obey g*(uj) = g{—uo). g' = g'(u>) 
is an even and g" = g"(uj) an odd function. This then implies 



9 



_ pr r- iu . (A 7) 

J 7T U l - UJ 2 



: ' analytic and free of singularities 
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A CAUSAL FUNCTIONS 



Now, only positive frequencies need to be considered. The counterpart to (A. 6) 



is 



Let us study an example. 
The Fourier transform of 

f(t) = 9(t)e- Tt (A.9) 

is 

sM = j^- . (A.10) 

Indeed, this Fourier transform has a singularity at uj = —iT in the lower half 
plane (we assume V > 0), and it respects g(u>) = —g*(—uj). The real and 
imaginary parts read 

= f4 • 9 " M = f4? ■ (A - U) 

they are even and odd functions of u>, respectively. 

_ f du 1 1 1 , . _ 

J TT U + tT U — id r — ltd 

is easily established, and the real part and imaginary part of this identity are 
the dispersion relations (A. 5) and (A. 6). 

By the way, (A.9) is the Green's function for a relaxation process. 

p + Tp = F(t) (A.13) 
is solved by 

p(t)= [ dr e~ TT F(t - t) . (A.14) 
J o 
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B Crystal symmetry 

The electromagnetic interaction, which governs the field of solid state physics, 
is invariant with respect to translations in space and time, rotations, space 
inversion, time reversal, and boosts. Let us concentrate on a body at rest, hence 
the boosts need not be considered any more. We have already investigated the 
consequences of invariance with respect to time reversal, so let us disregard 
this aspect here as well. We remain with spatial translations, rotations, and 
inversion. We first discuss why a symmetric Hamiltonian may give rise to an 
unsymmetric equilibrium state. This phenomenon, the spontaneous breakdown 
of symmetry, is widespread, in particular, when a system is very large, like a 
carbon hydrogen molecule or an ideally infinite crystal. 

B.l Spontaneous symmetry breaking 

The Hamiltonian H of ordinary matter is invariant with respect to transla- 
tion, rotation, and inversion. Consequently, the equilibrium, or Gibbs state, 
which is a function of the Hamiltonian, should also be invariant with respect to 
translations, rotations, and inversion. This is obviously wrong. Crystals have 
preferred directions (crystal axes), and they need not have an inversion center. 
In general, the symmetry of a theory does not imply the symmetry of the 
ground state of a particular system. Symmetry may be broken spontaneously. 
Just think about the planetary system. Although there is rotational invariance, 
the planets do not move on spheres (how should they?) or circles. Rotational 
symmetry merely says that a planetary system which is tilted by a certain angle 
would be possible as well. 

However, quantum mechanics teaches something else. The hydrogen atom, for 
example, is not a minute planetary system. All possible electron orbits interfere 
in such a way that the ground state is truly spherically symmetric. 
Another, less trivial example is the ground state of the ammonia molecule 
NH3. The three protons form an equilateral triangle. The nitrogen ion is 
either above or below. In fact, the ground state is a symmetrical superposition 
of above/below, and its energy is lower, by a tiny amount, than the energy 
of the anti-symmetrical superposition. The energy difference foujQ defines the 
microwave frequency standard (fo = 23.87012 MHz). 

The larger the molecule, the less probable is a transition between a state and its 
mirror-image. The right-handed version of a glucose molecule may, in principle, 
become a left handed version. The true ground state is a symmetric combination 
of both, but since so many nuclei would have to change position simultaneously, 
the transition time exceeds the age of the universe by orders of magnitude. 
Hence, if produced by a right handed plant, a dextrose molecule will stay in a 
state which defies the principle of inversion symmetry. 

The same applies to crystals. For instance, if LiNbO"3 is grown, an electric 
current defines a preferred direction. Once the crystal cools down, it keeps its 
c-axis orientation. Not even the largest external electric field may revert the 
crystal's polarization. 
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B CRYSTAL SYMMETRY 



B.2 Symmetry groups 

One has to distinguish between translations and point transformations. A crys- 
tal lattice consists of translated copies of a unit cell. A point symmetry element 
sends each ion of the cell to another position such that the transformed cell looks 
the same. The point symmetry elements can be an inversion center, a mirror 
plane, a 1, 2, 3, 4, or 6-fold rotation axis, or a 1, 2, 3, 4, or 6-fold rotation axis 
with inversion. These point symmetry elements form groups, altogether 32, 
each describing a crystal symmetry class. There is an internationally accepted 
short hand notation system. For example, 3m says that there is a three-fold 
rotation axis, without inversion, and a mirror plane such that the rotation axis 
lies in that plane. 3m would indicate an additional reflection symmetry with 
respect to a plane orthogonal to the rotation axis. Consult the standard work 
by Nye 16 for details. 

B.3 A case study 

Let us discuss lithium niobate LiNb03, a non-centric crystal of class 3m. This 

symmetry class allows proper tensors of rank 3. 

The symmetry group consists of the identity 7, a reflection 
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and a 120° rotation 
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The entire group 3m is made up of {7, D~,R,R 1 , 77 R, 7777}, because of R 2 
7? _1 , 77 2 = 7 and 777? 1 = 7277. The multiplication table is 
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16 J. F. Nye, Physical Properties of Crystals, their Representation by Tensors and Matrices; 
Oxford University Press 



B.3 A case study 
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We denote by c the polar three-fold symmetry axis (a vector of unit length) 
and by it, v, w three unit vectors in a plane orthogonal to it, with angles of 
120° between them. We identify ii = x and choose y orthogonal to it. This 
means 

i, = -^l x + ^ly and w = -^x-^y . (B.4) 

x — > —x should be a symmetry as well as a permutation of ii, v, and w. c — ► — c 
must not be a symmetry. After all, we deal with 3m, not 3m. We look for a 
third rank tensor d^k = dikj which would describe the Pockels effect or second 
harmonic generation. 

= ClCj c k (B.5) 
fulfills all requirements. 

Ujiih + VjVk + WjWk is a symmetric second rank tensor with three-fold rotation 
symmetry. A short calculation shows that it is proportional to XiXj + yiijj- 
Consequently, 

D ijl = Ci^jXk + VjVk) (B.6) 
is another admissible rank 3 tensor. A third possibility is 

(3) 

D ijk = ^i(cjX k + XjC k ) + yi(cjyk + VjC k ) . (B.7) 
Another tensor, namely UiUjiik + ViijjVk + WiWjWk, proportional to 

Xi(xjX k - yjyk) - yiixjyk + yjXk) , (B.8) 



is not acceptable, because it is antisymmetric with respect to x — > —x. How- 
ever, we only need to replace u by c x ii, and the same for v and w which 
amounts to interchanging x and y. Instead of (B.8) we now obtain the fourth 
tensor 

D ijl = Viivm - XjX k ) - Xi(xjy k + VjXk) . (B.9) 

Not that the first three tensors are invariant with respect to x <-> y, so c x u 
instead of u etc. does not produce new tensors. For the point group 3m there 
are just four linearly independent tensors of rank 3. 
The most general rank 3 tensor with 3m symmetry may we written as 

4 

d ijk = J2d r D^ k (B.10) 

r=l 

with four different invariants d r which depend on the effect and the material 
under study. 
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B CRYSTAL SYMMETRY 



We list the non- vanishing tensor elements: 



di 
d 2 
d 3 
d^ 



G?333 , 

dsn = d%22 , 

dl31 = ^113 = ^223 = ^232 , 
^222 = —d2U = — d\\2 = —d\21 



(B.ll) 
(B.12) 
(B.13) 
(B.14) 



A word of caution. Some authors represent a symmetric matrix by a six com- 
ponent vector: 



T ■ — 
j-ij — 



' T\ T 6 T 5 N 

T 6 T 2 T 4 
v T 5 T 4 T 3 j 



(B.15) 



When summing over indices, they count the off-diagonal index pairs, such as 
6 = (1, 2) = (2, 1), only once. Consequently you will find factors 2 in some such 
tables. The rank 3 tensor dijk is represented by a 3 x 6 matrix du where i runs 
from 1 to 3 and t from 1 to 6. In our case, (B.14) might read d\§ = — 2^22 • 
Nye 1 follows this convention. 



17 Physical properties of Crystals, loc. cit. 
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C Glossary 

Birefringence In the absence of gyrotropy the refractive part of the suscepti- 
bility tensor Xij , and with it the permittivity tensor = 5ij + Xij are real and 
symmetric. There is a Cartesian coordinate system such that the permittivity is 
diagonal, = <5y. A plane wave, linearly polarized along the \ th coordinate 
axis, propagates with main refractive index n = V eW. If all three eigenvalues 
of the permittivity tensor coincide, the medium is optically isotropic. If two 
are equal and differ from the third, there is a preferred axis. Light propagating 
along this optical axis may be polarized arbitrarily, but perpendicular to the 
optical axis. The corresponding main refractive index is called ordinary. Light 
propagating perpendicular to the optical axis is characterized by the main ex- 
traordinary refractive index. Crystals with an optical axis, or uniaxial crystals, 
are birefringent. If all three eigenvalues of the permittivity tensor differ, one 
speaks of a biaxial crystal. 

Cotton- Mouton effect A quasi-static external magnetic field, or the satu- 
ration magnetization of a ferro- resp. ferri-magnetic crystal, also contributes 
in second order to the susceptibility (Cotton-Mouton effect). The suscepti- 
bility change Ax^j 11 = x^ k( M^Mi is not bound to special crystal properties, 
Xifki fiijfikt is allowed. Because it is quadratic in the magnetization, the 
Cotton-Mouton effect is reciprocal, in contrast to the Faraday effect. 

Dissipation-fluctuation theorem If a system is disturbed by a time depen- 
dent addition -X(t)V to the Hamiltonian, the absorption of energy is described 
by the imaginary part of the corresponding susceptibility. Fluctuations of the 
equilibrium state are measured by the spectral density. Both functions, the 
dissipative part of the susceptibility and the spectral density, are intimately 
related, because of the Kubo-Martin-Schwinger identity. In particular, the 
Wiener- Khintchin theorem allows the proof of a weak form of the second law 
of thermodynamics. 

Faraday effect An external quasi-static magnetic field B causes a change of 
the susceptibilty by &X% = i> e ijk9ki where the gyration vector g (see gyrotropy) 
depends linearly on the magnetic field strength, or the saturation magnetiza- 
tion. If linearly polarized light passes a medium with Faraday effect, its po- 
larization vector is rotated by an angle 0L, where L is the path length within 
the medium and the specific Faraday rotation constant. If the light beam is 
reflected and passes the medium a second time, but backward, the polarization 
vector is rotated by the same amount, such that both angles add. This property 
allows to build an optical isolator. 

Gibbs state The equilibrium states of matter G(T, A) oc exp —H(\)/k-B,T, or 
Gibbs states, depend on temperature T and on external parameters A appearing 
in the Hamiltonian H. Besides chemical composition these may be mechanical 
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stress, an external electric field, or an external magnetic field. If the exter- 
nal parameters or the temperature of the system's environment change slowly 
enough, then Wt ~ G(T t , At) is a good approximation. The state Wt is always 
very close to an equilibrium state, and we speak of a reversible process. If the 
external parameters vary rapidly, the system is always away from equilibrium, 
and we speak of an irreversible process. 

Gyrotropy A gyrotropic contribution to the susceptibility is formally de- 
scribed by AxfJ = ieijkdk- 9 1S the gyrotropy vector. If light propagates in 
the direction of g, the eigenmodes are circularly polarized waves with different 
refractive indices. When a linearly polarized wave enters a gyroscopic medium, 
it is split into circularly polarized waves which propagate with different phase 
velocities. If this beam leaves the medium, the circularly polarized waves are 
recombined into a linearly polarized wave. The polarization vector of this out- 
going wave is rotated by an angle @L, where L is the path length within the 
gyroscopic medium. Note that the sign of rotation, right (clockwise) or left 
(anti-clockwise), is judged by an observer facing the light source. Optical ac- 
tivity and the Faraday effect cause gyrotropy. 

Hall effect If charged particles within a solid move in a static magnetic field, 
the current contains a component which is proportional to the driving elec- 
tric field, proportional to the quasi-static induction field, and perpendicular 
to both. These properties are a consequence of Onsager's theorem for kinetic 
coefficients, such as the tensor of electric conductivity. The sign of the Hall 
constant depends on whether electrons or electron wholes dominate the charge 
transport mechanism. 

Interaction picture The state W is defined by preparing the system un- 
der study in a well-described manner. An observable M is a class of equiva- 
lent measuring procedures which, for all states, show the same results. Both 
are represented, in conventional quantum theory, by self-adjoint linear oper- 
ators mapping a Hilbert space into itself. Time enters the game as follows. 
After preparing the state W and before measuring M, one may wait for a 
time t. This either defines a state Wt (Schrodinger picture) or an observable 
Mt (Heisenberg picture). Waiting is described by the unitary waiting operator 
Ut = exp (—itH/h), an exponential function of time t and an observable H, the 
energy, or Hamiltonian. 

The interaction picture is in-between. Often the Hamiltonian is the sum of 
a manageable part H and a perturbation V. It is advisable to resort to the 
Heisenberg picture with H. Then the Schrodinger equation for states is driven 
by V only, giving rise to a power series expansion in V. 

Kerr effect The second order electro-optic Kerr effect Ax^ e = xjjke^kEe is 
allowed even for otherwise isotropic media. It is usually much weaker than the 
Pockels effect, if present. 
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Kramers- Kronig relation The Hermitian part x[j °f the susceptibility ten- 
sor describes refraction, the anti-Hermitian contribution ix"j absorption. The 
refractive part at a certain angular frequency oo is an integral over all frequencies 
u of the absorptive part, weighted by l/(ui — u). There is no refraction without 
absorption, although possibly in another frequency domain. And: refractive 
indices depend necessarily on angular frequency. The Kramers-Kronig relation 
is a consequence retarded response functions. 

Kubo-Martin-Schwinger identity Formally, the Gibbs state is equivalent 
to the time translation operator with the inverse temperature as an imaginary 
time. This results in (BA(z)) = (A(z — ih(3)B) for an arbitrary imaginary 
time z. A and B are observables, A(t) is the time-translated observable, and 
A(z) a proper analytic continuation. The KMS identity is an important step 
in deriving the dissipation-fluctuation theorem. 

Linear response A system, originally in an equilibrium, or Gibbs state, may 
be perturbed by an explicitly time-dependent addition — A(t) V to the otherwise 
constant Hamiltonian H. The equation of motion and the initial condition may 
be formulated as an integral equation which can be expanded into a power 
series in A. The first non-trivial contribution describes the response of the 
system to such perturbations which is linear in the driving force A = A(t). 
By construction, the linear response, as described by response functions, is 
retarded. 

Onsager relations A symmetry relation for the generalized susceptibility 
tensor based on the invariance of physical laws with respect to time reversal. 
The susceptibility tensor is symmetric provided the direction of an external 
magnetic field is reversed. It can be extended to the tensor of kinetic coefficients 
describing the linear relationship between generalized fluxes and driving forces. 
Onsager relations are relevant for the tensorial properties of various effects, such 
as birefringence, the Pockels effect, Kerr effect, Faraday effect, Cotton-Mouton 
effect, optical activity, and the Hall effect. 

Optical activity The susceptibility depends on angular frequency and on the 
wave vector, Xij = Xijfa, q). The dependence on q is weak, and Xij(v) suffices 
for most applications. In some cases, however, the first order of an expansion 
in q will result in an detectable effect, namely in a rotation of the polarization 
vector along the propagation path. This gyrotropic effect is known as optical 
activity. It will occur only if the material is left-right hand asymmetric, such 
as quartz or a solution of dextrose. 

Optical axis A direction such that a beam propagating along it has a re- 
fractive index which does not depend on polarization. For optically isotropic 
media, any direction is an optical axis. If two eigenvalues of the permittivity 
tensor coincide, but differ from the third, there is but one optical axis. If all 
three eigenvalues are different, one has two optical axes. 
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Optical isolator A bulk optical isolator consists of a polarizer, a crystal with 
Faraday effect rotating the polarization by 45°, and another polarizer being 
rotated by 45° with respect to the first. If light which has passed the isolator 
is reflected and possibly de-polarized, the polarization of the reflected wave 
passing the second polarizer is rotated by another 45°; it is therefore blocked 
by the first polarizer. Optical isolators are required to protect a laser from its 
own light. An integrated optical isolator is the goal of intensive research effort. 

Pockels effect A linear opto-electric effect Ax^ c = xljk^k is possible for 
ferro-electric or piezo-electric crystals which possess no inversion symmetry cen- 
ter. This Pockels effect allows to switch or modulate light beams. It is usually 
much stronger than the Kerr effect. 

Refractive index Plane harmonic waves are characterized by an angular 
frequency u> and a wave vector q. With w (a unit vector) as the direction 
of propagation and ko = uj/c one may write q = nk^w. This defines the 
refractive index n of the wave. If the susceptibilty, consequently the permittivity 
tensor = Sij + Xij 1S rea l an d symmetric (no absorption, no gyrotropy), one 
may choose a Cartesian coordinate system such that e becomes diagonal. The 
square roots of the diagonal entries (eigenvalues) are sometimes called main 
refractive indices. They are refractive indices if the wave is polarized along the 
corresponding main axis. 

Response function If H t = H — \{t)V drives a process, than the expecta- 
tion value of an observable M at time t is the equilibrium expectation value 
plus a retarded integral over the driving force A = X(t). The Green's, or re- 
sponse function r(r) is proportional to the equilibrium expectation value of 
the commutator [M(r),V]. The Fourier transform over positive ages r is a 
susceptibility. 

Second harmonic generation The quadratic response of a medium to per- 
turbations of angular frequency to produces a polarization of double frequency. 
Only crystals without inversion center show this effect, such as lithium niobate. 
Efficient frequency doubling requires phase matching. 

Susceptibility If the Gibbs, or equilibrium state is disturbed by a time de- 
pendent addition —\{t)V to the Hamiltonian, the expectation values of ob- 
servables M vary with time, at lowest order with the same frequency as the 
perturbation. The Fourier transform of (M)t — (M) is linearly related with 
the Fourier transform A, the susceptibility x = x( u ) being the constant of pro- 
portionality. Susceptibilities are Fourier transforms of corresponding response 
functions r(r) over positive ages r. 

Wiener-Khintchin theorem The process t — ► M{t) in a stationary state 
is described by the time auto-correlation function K{t) indicating how much 
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M(t) is still correlated with M(0). The Fourier transform of K{t), the spectral 
density S(uj), is never negative. A very short correlation time results in an 
almost constant spectral density which qualifies the process t — ► M{t) as white 
noise. The Wiener-Khintchin and the dissipation fluctuation theorem allow to 
prove the second law of thermodynamics. 



